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Abstract. We consider two-dimensional degenerate electron gas in the presence of 
perpendicular random magnetic field. The magnetic field disorder which is assumed 
to be gaussian is characterized by two parameters. The first is proportional to the 
amplitude of local magnetic field fluctuations. The second one characterizes the 
disorder on the longer scale - we call it the screening length. Using Kubo formula 
for the conductivity we have found a class of diagrams which leads to the long time 
tails in the current - current correlation function. For short times comparing to 
the difi^usion time corresponding to the screening length this function behaves like 
logarithm of time, for longer times it decays like what of course puts in question 
the diffusive character of the behavior of the charged particle in random magnetic 
field widely assumed in the literature. 



PACS numbers: 05.60 -fw, 72.10 -d 

The theoretical description of two-dimensional, quantum, noninteracting, 
degenerate electron gas moving in a spatially random magnetic field attracts nowadays 
considerable attention in view of its possible relation to theory of Quantum Hall 
Effect near the filling factor ly — 1/2 proposed by Halperin and collaborators [Q. The 
characteristic feature of these systems behavior is their temporal and spatial dispersion 
reflecting underlying microscopic motion nonlocality ||^. Progress in developing 
general theory capable of describing this dispersion in degenerate electron gas is rather 
slow. Attempts were directed, therefore, to analysis of relatively simpler problem, 
namely the particle diffusion. In the recent work Aronov, Mirlin and Wolfle |^ 
proposed the derivation of the zero frequency diffusion coefficient for a two-dimensional 
degenerate electron gas in random magnetic field. In the present paper we propose 
a next step in the direction of the full analysis by using the Kubo formula for the 
conductivity. Analyzing a particular class of diagrams we have found the long time 
tails in the current current correlation function. 

We consider a gas of noninteracting, charged, spinless fermions moving on a 
two-dimensional plane in the presence of a spatially random magnetic field B(r) 
perpendicular to that plane. The single particle Hamiltonian reads: 

H=±(^W-'-Ay (1) 
2m \ I c I 



where A is the magnetic field vector potential. 
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The paramount difficulty is now the choice of statistical ensemble for the random 
magnetic field. The seemingly obvious, and natural, choice of Gaussian distributed 
magnetic field with the second moment {Bz{y)Bz{v')) ~ 5{y — r') results in infrared 
divergencies in the lowest order terms in the perturbation theory This is due 
to the long range vector potential correlations reflecting the fact that above Gaussian 
distribution for magnetic field permits the field configurations in which arbitrary large 
magnetic flux is penetrating the plane. The mean square value of magnetic flux inside 
a region is proportional to its area, and this invalidates perturbational approach to 
the problem. 

To avoid this we consider a model of the magnetic field disorder in which the above 
mentioned configurations are explicitly suppressed. To this effect we assume that the 
components of the vector potential Ai{v) (i = x,y) are the Gaussian distributed 
quantities with zero mean, and which second order cumulant, in the Fourier space, is 
given as 

{A,{q)A,iq')) = (2^)2<5(q + q') (^^^\-^Sf^{q) (2) 

where ^^(q) = Sij — qiqj/o^, and 7 and are some constants. This ensemble is 
closely related to that used by Bausch, Schmitz and Turski Q in their analysis of the 
classical particle diffusion in crystals with randomly distributed topological defects. It 
is also identical to the one proposed by Aronov, Mirlin and Wolfle Q . In their paper 
7 was proportional to the fluctuations of the magnetic field and ^ was introduced 
at the intermediate steps of calculations as a "soft cut-off' in order to avoid the 
above mentioned infrared divergencies in the perturbation theory. They didn't discuss 
the meaning of this parameter, because it drops out from their final formula for the 
diffusion coefficient. This is true for the calculations in the lowest order in 7, as we 
shall see the higher order terms contain /i in the very nontrivial way. 

One can attach a clear physical meaning to the parameters 7 and /i in Eq. 
by analyzing an analogy between the above discussed model and another model of 
randomness in which these coefficients appear in a quite natural way. This auxiliary 
model has a well defined physical interpretation and long and short wavelength 
behaviors of both these models are identical. 

Imagine two identical, infinitely thin solenoids (Aharonov - Bohm flux tubes) 
piercing our plane at points separated by a distance £. The magnetic flux in one 
solenoid is +$ and in the second — <&. For sake of definiteness we shall call such an 
object a flux dipole. On the two-dimensional plane the dipole is characterized by its 
position, i. e. the coordinates of its "center of mass" and the angle tp between the 
vector joining the flux tubes and, say, the x axis. 

Now we assume that we have many such dipoles penetrating the plane and let 
the average dipole density be n. The position of the i'th dipole Ri may be arbitrary 
and the angle Lpi is homogeneously distributed in the interval (0, 27r). Notice that 
by construction our disorder ensemble allows only for configurations in which the 
total magnetic fiux through the plane equals zero. For such kind of the magnetic 
disorder we obtain the following expression for the second moments of the Fourier 
space components of the vector potential field Ai (q) : 

(A,(q)A,(q')> - C(q,q') = (271)2^ (q + q')2n (^^-^ g(|q|)5^^(q) (3) 

where 5(|q|) — (1 — Jo(|q|^))/q^, Jo{z) is the zero order Bessel function and <f>o = e/fic. 
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The ensemble average denotes here 



i=l 







where N s the number of dipoles and S is the area of the system {n = N/S). 

We shah also assume that the magnetic flux $ <C <&o, thus we may neglect higher 
order cumulants which are of the order n<I>'' with (r > 4) . Alternatively we may assume 
that the fluxes in different dipoles are mutually independent, gaussially distributed 
quantities with zero mean and the second moment equal to $^ . In this case the higher 
order cumulants also vanish. 

Now discuss some properties of the introduced disorder. If the condition y/n£ ^ 1 
is satisfled then on the length scales shorter than £ (we call £ the screening length in 
the following) the magnetic fleld is spatially uncorrelated because the magnetic fluxes 
in randomly placed solenoids are mutually independent. This follows also from the 
behavior of C(q, q') for q£ ^ 1. In this limit 5(|q|) ~^ 1/q^, as for the delta correlated 
magnetic field. In the opposite limit, \q\£ — > , the function gdq]) const and we 
have the case of spatially uncorrelated vector potential field. Comparing the short and 
long wavelength behavior of the models defined by Eq. (||) and Eq. (||) we see that 
they are indeed identical provided 7 = 2n($/<I>o)^ and /i = 2/£. Whichever of the 
limiting behavior of the vector potential correlation is pertinent to an actual physical 
system depends on the value of the Fermi wavelength for the electron gas. In this 
paper we are interested in situation when the Fermi wavelength Xp = 2'K/kp is much 
smaller than £, i.e. kp£ ^ 1. 

In the zero temperature and to limit Kubo formula for the real part of the 
conductivity reads: 

Rea^,{iu, k = 0) = ^^^M Tr [V:,Gr{Ep + u)V^ [Ga{Ep) - Gr(Ep] 

where Vx — \dx ~ ^A^iv) is proportional to the x'th component of the velocity 
operator and Gria{E') are retarded and advanced propagators corresponding to the 
one particle hamiltonian (|^) . In position representation they have the form: 

Gn/A{E,r,r)-}_^^—-^ (6) 

where <y3„(r) and e„ are eigenfunctions and corresponding eigenenergies. In the equa- 
tion (||) Tr denotes integration over r and S is the system's area. 

In our analysis we apply the standard method of the averaging over the disorder 
lH, I?). Every Green's function in Eq. is expanded in the power series of the 
components of the vector potential which we treat as a perturbation. In the second 
step every term of the resulting series is averaged according to the assumed properties 
of the magnetic field disorder. The exact summation of that series is impossible. What 
is usually done is the extraction and summation of certain classes of diagrams which 
are believed to give the most important contribution to the final result. The criterium 
of "importance" is based on the behavior (an index of divergence) of a corresponding 
analytical expression in the limit w ^ In the present case the most divergent 
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diagrams may be divided in two classes. The first class leads to renormalization 
of the one particle Green's function (the averaged Green's function). The second 
class consists of the ladder diagrams. These two classes when summed up give the so 
called Boltzmann value of the conductivity , or equivalently the Boltzmann diffusion 
constant, if one uses the density-density response function instead of Eq. and 
extracts its low frequency and long wavelength behavior ^ . 

This procedure well established and widely used in electrostatic potential disorder 
case, however, suffers here, from certain difficulty. Namely, the self energy in the 
averaged Green's function calculated for the disorder ensemble defined by Eq. (^) is 
infinite in the limit /i = j ^0. The averaged Green's function in the limit kl >> 1 
equals 

(G^/^ {E, k))^l/{E- fc V2m ± iujo/2) (7) 

where ujq = / ^q)^ kl / {2m) and even for very small density of dipoles n or for 
very small amplitude of magnetic field fluctuations 7 the self energy operator becomes 
dominant for large I. Of course we realize that the self energy in the averaged Green's 
function has no direct physical meaning (due to the fact that the Greens function is 
not a gauge invariant object) thus this infinity may be of no importance. Such a point 
of view seems to be justified, because in the final expression for the diffusion constant 
(a gauge invariant quantity) calculated in the ladder diagram approximation the limit 

— > may be performed and the diffusion coefficient is finite. 
On the other hand the methods of calculations of the diffusion coefficient (for example 
by solving the Dyson equation in the ladder diagram approximation [^) depend on the 
assumption that loq ^ Ep- Thus the question arises whether the discussed difficulties 
concerning the breaking up the perturbation theory are spurious and arising on the 
intermediate steps of calculations only do not influence the final results or they are of 
physical origin and one should be very careful in dealing with them. 
We think that the simplest method to answer this question is to perform the analysis 
of the diagrams' divergencies similar to that presented by Kirkpatrick and Dorfman 
for the electrostatic potential disorder problem, i. e. to analyze the divergencies 
with respect to co in the given order of the expansion parameter 7. It turns out that 
it is enough to consider diagrams of the first order to notice some very interesting 
features. These diagrams are shown in Fig. 1. Here the continuous lines correspond to 
the unperturbed Green's functions G^{E, k) = 1 / [E — k"^ / 2m±iri) and the broken line 
represents the averaging. The number of diagrams is much larger comparing to the 
case of electrostatic potential fluctuations case. This is due to first: the appearance of 
the second power of the vector potential in the hamiltonian - this term must be kept in 
order to conserve the gauge invariance, second: the appearance of the gauge invariant 
derivatives V in Eq. (^). The origin of the first five diagrams la - le is the expansion 
of the Green's function in powers of the vector potential. In particular, diagrams Id 
and le correspond to the term in the hamiltonian. The diagrams If - li comes 
from the expansion of the Green's function and from covariant derivatives, and the 
last diagram Ij represents the contraction of two A's in two covariant derivatives in 
Eq. (I). 

The first three diagrams are the lowest order diagrams belonging to the above discussed 
classes and although each of them separately is divergent in the limit / — > cx), the sum 
of them has well defined limit and it constitutes the first order term in the density 
expansion for the frequency dependent conductivity. That sum behaves like 7/0;^. 
The diagrams Id and le do not give contribution to the real part of the conductivity. 
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The contributions from the next four diagrams cancel. The last diagram Ij has an 
interesting behavior. Its analytical expression is proportional to the integral 

fc2 + - 2fcgcos(6') + ^2 > 



where k = ^2m(E + fiLo^jfi and q — \/2raE/fi. We see that for /i 7^ it is finite 
for w = 0. However for /i = it behaves like 7/w. Such a behavior suggests possible 
relationship between uj and ^ when both of them are small. 

Before proceeding further let us notice that up to now all the calculations have been 
performed in the Coulomb, symmetric gauge. The same result however is obtained in 
the Landau gauge = 0. In that case all the diagrams originating from the covariant 
derivatives vanish and we are left with first five diagrams la-le. Again the diagrams 
Id - le do not give contribution and for /i = the ^ /u) divergence obtained previously 
from diagram Ij is contained now in the sum of first three diagrams la - Ic. Of course 
this is nothing special, it should be so because the expression (H) is gauge invariant. 

The natural generalization of the diagram Ij is presented at Fig. 2. Contrary to 
Fig. 1 the solid lines now represent the dressed Green's functions. It is interesting that 
assuming the diffusive behavior of the electron gas we may calculate the sum of these 
two diagrams exactly. Namely, without the broken lines the sum of two diagrams 
is closely related to the density - density correlation function, more precisely it is 
proportional to the Kubo relaxation function which for q — > and w — > behaves 
like 

$(q,c.) = -AfII^ (9) 

That is why the analytical expression corresponding to the sum of the diagrams at 
Fig. 2 may be calculated and the final result reads: 



Dp 



The Eq. ^ has no limit for lu ^ {ReSa ln(u;)). 

It is interesting to write the sum of diagrams from Fig. 2 in the " time representation" 



6a = i— / dte"^'F(t) (11) 

IGn^m Jo 

where 

/•CO 

F{t) = / dq^^e-^'^'' (12) 

is proportional to the current - current correlation function. From the above we see 
that for small times comparing to the diffusion time corresponding to the screening 
length {Dfi'^t ^ 1) F{t) ~ ln{Dfj,'^t) and for asymptotically large times {Dfi'^t 1) 

The obtained result is interesting and not entirely clear. First, in conjunction with 
the diagrammatic analysis presented earlier it suggests that the screening length or the 
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" soft cut-ofF in the terminology of Aronov, Mirlin and Wolfle has important physical 
meaning determining time characteristics for the system. The power law decay of the 
function F(t) is of the great interest. The behavior suggests that the long time 
tails are of quantum origin j^], for the classical Lorentz gas model would give 
behavior. To fully appreciate whether this decay is indeed of a quantum mechanical 
origin requires a detailed analysis of the velocity - velocity correlation function for the 
classical particle in random magnetic field. To the best of our best knowledge this is 
not available at the moment. 

We should now return briefly to the discussion of the applicability of the pertur- 
bation theory. According to Eq. perturbation theory is applicable if hujQ <C Ep = 
"flXp <SC 1- Consider now a part of the system, linear dimension of which is i? ~ j-'^/^ 
(i? <C 0- The typical magnetic field flux contained in that part of the system $ ^ $o, 
consequently the typical magnetic field B ^ <&o/i?^ ~ 7$o- For a particle of energy 
Ef the corresponding cyclotron radius n = "^^J^^" - - I/7AF. Thus the 

condition ^IXp ^ 1 may be rewritten as 3> ^ what means that the cyclotron radius 
corresponding to the local magnetic field fluctuation should be much larger then the 
screening length. 

In conclusion we have studied the behavior of the two dimensional degenerate 
electron gas in the presence of the perpendicular random magnetic field. Analyzing 
the lowest order diagrams we have pointed out the importance of the proper choice 
of the magnetic field disorder ensemble. Secondly, taking into account certain class 
of diagrams we have obtained the long time tails in the current - current correlation 
function. 
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FIGURE CAPTIONS 

FIG. 1. Diagrams for the first order perturbation theory for the conductivity. The 
continuous lines correspond to the nonperturbed Green's functions, the broken lines 
to the averaging over the disorder, respectively. 

FIG. 2. Diagrams leading to the nonanalytical corrections to the conductivity. Here 
the solid lines represent the dressed (averaged) Green's functions. 



